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Abstract 

We consider a BCS-type model in the spin formalism and argue that the struc- 
ture of the interaction provides a mechanism for control over directions of the spin 
S other than S z , which is being controlled via the conventional chemical potential. 
We also find the conditions for the appearance of a high-T c superconducting phase. 
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Introduction 

Twenty years after the discovery of high-temperature superconductivity [1] there is still 
neither consensus nor clear understanding of the mechanism or mechanisms which are 
behind this exciting and with innumerable practical applications phenomenon. The initial 
discussions (see, e.g. [2]) have led to the formation of some main conceptual stream, as 
presented in [3], however other view points are continuously being argued, just to mention 
a recent one [4] . 

In [5] we proposed a combination of a BCS and a mean-field Hamiltonian where the 
transition temperature could become arbitrarily high. This happened without increasing 
the interaction indefinitely but by a small denominator. In this note we investigate this 
effect more closely and find that another important ingredient is a chemical potential 
which breaks the electron conservation. We give a model for this phenomenon by the in- 
teraction with a reservoir of quasi-particles which do not have a definite electron number. 
Since such objects play an important role in the theory of the Josephson currents [6], we 
think that this possibility is not purely academic. 

To avoid a terminological misunderstanding, we recall that the (quantum-mechanical) 
mean-field theory and the BCS-theory of superconductivity correspond to essentially dif- 
ferent physical situations. A mean-field theory means that the particle density p(x) = 
ip*(x)ip(x) (in second quantization) tends to a c-number in a suitable scaling limit. With 
an appropriate smearing, from the operator-valued distribution p(x) an unbounded oper- 
ator is being produced, so that the best to be strived for remains the strong resolvent con- 
vergence in a representation where the macroscopic density is built in. In the BCS-theory 
pairs of creation operators with opposite momentum ip*{k) ip*(—k) tend to c-numbers, so 
the correlations required in both cases, seem to be quite different. The main result in [5] 
was that both types of correlations may well co-exist in certain regions of the parameter 
space (temperature, chemical potential, relative values of the two coupling constants) 
and this appears to be the case in the KMS-state of the equivalent approximating (Bo- 
goliubov) Hamiltonian H B , two Hamiltonians being considered as equivalent if they lead 
to one and the same time evolution of the local observables [7]. 

In what follows, we generalize the original BCS model in the most natural way, namely 
by augmenting it with the missing mean-field interaction components. We show that this 
provides a mechanism for control over directions of the spin ~3 other than S z , which is 
being controlled via the conventional chemical potential. 

1 The degenerate BCS Hamiltonian 

The initial quartic BCS Hamiltonian is mainly known in terms of fermionic creation and 
annihilation operators [8] 

H = Y,(^k - p){a\ ik a 1tk + a\ tk a ltk ) + ^ V k , k ,a\ k ,a\_ k ,a^ - k a ltk . (1.1) 

h h h> 

fx, i\j ■ (V 
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It involves however only the algebra generated by the pair-operators a^-kd^k (observe 
ajai + aja^ = [ajaj, aja^ + 1), so we shall only be concerned with them and shall represent 

them by spin matrices a \j a \~j —> < J j+ = (&j,x + i a j,y) l^i j — 1, . . . , N. As a weak 
interaction can only scratch the Fermi surface, we take oju = us, \fk and incorporate the 
latter into ji. Finally, we set Vk,w = —2X B /N Vk,k'. With the notation ~3 = at, 
Hamiltonian (1.1) becomes equivalent to 

H = -^{Sl + Sl)-^S z . (1.2) 

In this form the Hamiltonian can be exactly diagonalized and the following steps are 

mathematically rigorous in the limit. 

We assume that the thermal state (A) = Tr Ae _/3H /Tr e _/3H is such that the length of 
is much bigger than the fluctuations around it, ((S — ( S^)) 2 ). This means that in the 

identity 

S 2 = (S- (S)) 2 + 2{S)S- (S) 2 

the first term is small compared to the second one. Since the last term is a c-number, 
Hamiltonian (1.2) becomes equivalent to the following one, linear in 

H B = -X B + ^S tf ) - »S M . (1.3) 

Of course, the original Hamiltonian is invariant under rotations around the z-axis, but the 
spin-vector S will point into some direction (with —fj,S z contribution taken into account) 
and we shall call this resulting spin direction Sb, that is Eq. (1.3) can be rewritten as 

H B = ~*b^S x - fiS z =: -W S B , (1.4) 

where 

W = ^/(\ B (S x )/N) 2 + ^ 2 >\fi\ (1.5) 

S B = bS x + VT=WS x , (1.6) 

and the mixing parameter of the Bogoliubov transformation [9] is defined through 

b = X B (S X )/NW . (1.7) 

This rotation can be inverted and if S± is in the x — z plane orthogonal to S B , we 
have (Figure 1) 

S x = bS B — Vl — b 2 S±. 

Since H = — W S B is the sum of N spins in the 5-direction, Eq.(1.4), the thermal 
expectation values are the usual ones 

%l = tanh^, (S ± ) = 0. (1.8) 
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The self-consistency of the system is expressed by the so-called "gap-equation" 

N 2T W N 2T K J 

This gap-equation has two solutions 

(A) a normal state 

(S x ) = VT (1.10) 

(B) a superconducting state 

Wn Wn 

(S x )^0, ^ = tanh^, (1.11) 
^T = ^ F m<T c , (1.12) 



2 V2T 

where the characteristic function < F(ct) < 1 is given by 

1 - a 2 /3, a -> 



. tanha y 

F(a) = 

a \ 



1/a, a — > oo 



(see Figure 2) and the critical value which the temperature T for no values of /x and A 
can exceed is 

T c <^. (1.13) 
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Figure 2: The characteristic function F(a): high-T region corresponds to small a. 



On Figure 3, the pure BSC-situation is shown: the plot of both sides of Eq.(l.ll), for 
T = A s /4, A B /2,3A B /4. The limit (1.13) becomes obvious. 

Relations (1.5), (1.12) in fact suggest some possibilities for high-T c generation. A 
realistic mechanism should result in a deviation from the proportionality relation (1.12) 
as far as the characteristic-function part is considered (as F(a) is always less than 1). 
Also, it should aim a modification of the quasiparticle dispersion relation (1.5), as is e.g. 
the case of the gossamer superconductor [10]. We rather target the appearance of an 
effective chemical potential, whose variation would provide a means of influence on the 
transition temperature. 



2 On the role of the chemical potential 

The chemical potential is a control parameter which adjusts the number of Cooper pairs, 
in our formalism S z . We shall now argue that the BCS-interaction gives us a handle to 
control also other directions of ~3 and in this way to clone the Josephson phase of R. 

Suppose our system interacts strongly with a superconducting reservoir R such that 
Eq.(1.3) holds for the ensemble: 



H = -T\ iS+S } N + N R + ) N + N R 



The cross-term 



(2.1) 



(2.2) 



N + N R 

induces control parameters which just copy on the system the situation in R, if R is 
dominant [7, 11, 12]. By the coupling with the reservoir, the !> -direction is cloned. We 
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shall call this spin-coaxial exchange. Thus, if the reservoir is in the normal state, is 
in the z-direction and we get the usual chemical potential. If R is superconducting, we get 



a chemical-potential coupled (a vS^'-term in the Hamiltonian) , which represents 
the quasi-particles — the elementary excitations of the superconductor. If the reservoir 
dictates a Josephson phase, e.g. in the x-direction, this means that in the cross-term (2.2) 
(S X R ^) ~ v. Correspondingly, if we treat the Hamiltonian with this additional term as 
before, it becomes 

H — uS x = -WS B , v > (2.3) 



with dispersion relation 



and mixing parameter 



The gap-equation reads: 



W 




N 



W 



W 2T 



(2.4) 



and the temperature at which the gap opens becomes 



1 / v_ 

1 = 2 \ (S X )/N 



\ W 



(2.5) 



So with this choice 

(i) expectedly, there is no normal phase, i.e. solution with (S x ) = 0; 
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(ii) VT there exists a solution with (S x ) > 0. 

Thus there is no phase transition W > 0. It is quenched, since the symmetry is broken 
externally and not spontaneously. 

With all this taken into account, for the description of a system which exhibits a phase 
transition towards a high-temperature superconducting phase, we study the following 
system: It has a BCS interaction of the form of (1.2) and a self-interaction ~ XmS z 
which can also be treated as a mean field for the states we are considering. Furthermore 
it is coupled to a particle reservoir which supplies a chemical potential /iS z . Finally, it 
is in interaction with a superconductor which imprints its phase by a term uSb — the 
component of S in suitable direction to be determined later. Thus we have a 4-dimensional 
parameter space V, two coupling constants - Xb and Am, and two chemical potentials. 
For each point in V there is a transition temperature T c (Xb, Am, A*, v) > and we want 
to exhibit a region 1Z e V, where this function is not bounded, that is to say that T c in 
this region can become arbitrarily high. 

In formulae we start with the Hamiltonian 



H 



2N 



which in the mean-field regime becomes 

tl = —ARJ X 



XmS; 



(S,) 



fxS z - vS B - 



(2.6) 



(2.7) 



N N 

Sb is the combination of S x and S z which prevails in the superconducting reservoir 
and should be a perfect match of the combination we are getting for our system. This 
sounds somewhat mysterious, as such a cloning seems to require some foresight from the 
external system. However, there is some redundance in the coefficient of S z which has 
contributions both from /i and v. We are supposed to be able to control \x and by giving 
part of it to v, we can use it to adjust the direction of Sb so that it coincide with the 
direction we shall obtain. Thus we write 



with 



H = (W + u) 



W 




Xm(Sz) 



(2.8) 



\ N 2 



Xm(Sz) 



'tiff 



N 2 



> Wff\ 



This Hamiltonian gives rise to the following system of coupled gap-equations: 

(Sx) 



N 
N 



NW 

Heff(Sz)/N 

W 



IT 



tanh 



W + u 
2T 



(2.9) 

(2.10) 
(2.11) 
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The system (2.10-11) has both solutions, corresponding to normal and to supercon- 
ducting phases, so with (S x ) = 0, resp. (S x ) ^ 0. In the latter case, from Eq.(2.10) the 
transition temperature at which the gap closes, (S x ) — > 0, is found to be 

Tc JM\^±M. (2 . 12) 

* |A*eff| 

As n e ff can be made arbitrarily small, this means that for given values of the parameters 
/x, v, \m and Xb, in a bounded region 1Z e V, the critical temperature can become 
arbitrarily high, as suggested in [5]. 

The second (coupled) gap-equation provides a relation between the model parameters 
that determines the relevant parameter range: 

(S z ) /i 



N Xb + Xm 

The existence of further order parameters is of severe importance for the physical 
content of the models under consideration [13]. Even in the simple model above, the 
presence of a second order parameter leads to an enrichment of the structure and to new 
effects. 



A B <0 X B >0 




Figure 4: Mean-field enhanced BCS: spin-coaxial exchange (plot of both sides 
of Eq.(2.15); the line thickness increases with T and A#). 



Let us discus the superconducting solution more in detail. In the mean-field enhanced 
model, for (S x ) ^ 0, Eq.(2.10) reduces to 

_ = tanh __ (2.13) 

We have chosen the positive eigenvalues of H, Eq.(1.5). This is not really a restriction, 
since the consideration of the opposite situation will give the conjugate picture. Thus, 



9 



tanh(...) and Xb must always have the same sign. Also, Eq.(2.9), the lower bound for 
the values of W is determined through the effective chemical potential in the ^-direction, 
/Xefr- Depending on the coupling of the system to the reservoir (the value and the sign of 
v), we are led to the following situations: 

(A) Spin-coaxial exchange, v > (Figure 4) 

• the BCS-coupling has to be attractive and stronger than the effective chemical 
potential, \b > \fJ>es\] 

• the solution (when existing) is uniquely determined; 

• as also seen from Eq.(2.13), the higher-temperature solutions require also stronger 
BCS coupling (the thickness of the lines increases with T, resp. with A#; the ad- 
missible solutions have to be to the right of the dashed line). 




Figure 5: Mean-field enhanced BCS: spin-anticoaxial exchange, (a) \fi c s\ < Mi 
(b) |/U cff | > \v\ (also here, the line thickness increases with T, resp. with \ B ; 
the admissible solutions have to be to the right of the dashed line) 



(B) Spin-anticoaxial exchange, v < 

In this case the relative values of \n c s\ and \u\ become of importance. 
When l/ieffl < \u\, (Figure 5a), 

• solutions with repulsive BCS-coupling are possible and uniquely defined; 

• \\b\, when \ B < 0, has to dominate the effective chemical potential, \/i e s\] 

• in the positive A^-coupling range, it can happen that the full system has none, one 
or two solutions. 
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When \(J, e ff\ > (Figure 5b), 

• only solutions with positive BCS-couplings are possible; 

• depending on the relations between the parameters — /z e fr, and v — encoded 
in the second gap-equation, the system can have none, one or two solutions. 

As Eq.(2.12) requires small values of |/i e fj| in order to achieve high transition temper- 
ature, this would correspond rather to the situations depicted on Figures 4 and 5a. 



3 Conclusions 

We considered a BCS-type model with two order parameters, whose solvability is encoded 
in two coupled gap equations and which exhibits a high-temperature superconducting 
phase. High-T c superconductivity models that are based on coupled gap equations, are 
known in the literature: such an approach is the one due to Eliashberg [14], see also [15] 
for recent analysis. There, the limitations on T c also disappear and are thus interpreted 
as artifacts of the Bogoliubov method. However we could not identify the underlying 
mechanism with the one described above. Some argumentation for a higher, compared 
to BCS, or unlimited transition temperature comes also from the line of considerations 
towards an unification of the BCS and BEC pictures [16]. Our model might be relevant 
here as well as it exhibits an off-diagonal long-range order, ODLRO [17, 18]. Recall that 
its existence is the basis for the Bose-Einstein condensation, however we are dealing here 
with a fermion system, so the model provides a framework for analysis of BEC in a Fermi 
gas [19, 20]. 
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